On the corresponding states law of the Yukawa fluid 
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We have analyzed the currently available simulation results, as well as performed some additional 
Monte Carlo simulation for the hard-core attractive Yukawa fluid in order to study its corresponding 
state behavior. We show that the values of reduced surface tension map onto the master curve, 
and a universal equation of state can be obtained in the wide range of the attractive Yukawa tail 
length after a certain re-scaling of the number density. Some comparisons with other nonconformal 
potentials are presented and discussed. 



PACS numbers: 



I. INTRODUCTION 



The investigation of colloid agglomeration, and phase 
stability has been a subject of long-standing theoreti- 
cal and practical interes^^^. It is well accepted that 
the range of the effective interaction between particles in 
different solvents determines the shape and location of 
the phase boundaries of colloidal dispersions^. From a 
theoretical point of view, the phase transition in colloid- 
solvent solution is analogous to that of the single compo- 
nent gas-liquid system. Thus, the most general approach 
is to assume a form for the interparticle effective pair po- 
tential. One of the most common selection is the Yukawa 
potential, 



u(r) 



o, if r < a, 

-e<7exp[—n(r — a)]/r, if a < r, 



(1) 



where parameters e and n define the depth and range 
of the potential, respectively; a is a particle diameter, 
which is usually taken as a unit length, a — 1 . 

The range of attraction in a pair potential cannot be 
changed for ordinary fluids, but it can be adjusted in 
colloid-solvent systems by adding, for instance, a nonad- 
sorbing polymers, ions, or other solutes^^. Therefore, 
upon varying the range of the interaction in the Yukawa 
potential (Q}, one can repoduce the behavior of some real 
systems. 

The phase diagram and interfacial properties of the 
hard-core attractive Yukawa (HAY) fluid have been in- 
vestigated by means of the theoretical tools of the 
liquid-state statistical mechanics^^^^, and computer 
simulation s 11 ! 12 ' 13 ! 14 ! 15 . The simulation for measuring 
any thermophysical property or phase equilibria for HAY 
fluid is, in most cases, time consuming and expensive due 
to long-range character of the potential. That is why 
HAY fluid with long attraction tails have been recently 
studied by applying an Ewald su m 15 ' 16 . 

On the other hand, it is well-known that in some cases 
the application of the corresponding states (CS) law is a 
very useful approach to avoid extensive computer calcula- 
tions as well as experimental measuring a 17 d 8 i 19 i 20 i 21 i 22 i 23 . 
According to the CS law the substances have unified 



equation of state 



Pr = F(T r , Pr ), 



(2) 



in reduced variables of pressure, Pr = P/P c , tempera- 
ture, Tr = T/T c , and number density, p^ = p/p c ] the 
subscript c denotes critical values, and T = ksT/e, as 
usual. In Eq. [2] F is a universal, but complicated, func- 
tion. However, it is well-known that the CS approach 
is valid only for conformal pair potentials, and Yukawa 
potential is not one of them, because its range of attrac- 
tion varies independently of the hard-core radius a. Tak- 
ing this into account, some researchers have applied the 
so-called extended CS law for a nonconformal potential, 
i.e. when Eq. [2] involves a third parameter—. Besides, 
there have been speculations based on experimental ob- 
servations and computer simulations for various noncon- 
formal model potentials that the gas-liquid coexistence 
curves may follow an extended CS theory, if they are 
plotted in terms of the reduced second virial coefficient, 
Bl = B 2 {T C )/B^ S (where B% s is the second virial co- 
efficient of hard sphere fluid), and the reduced density 
pp2&. Such considerations were based on the assumption 
that second virial coefficient remains practically constant 
at the critical point. 

However, more recently it has been shown, that B^ 
varies with the range of attractions, in contrast to the 
common belief that it remains practically constant for 
different pair potentials*^. The results of B^ for the 
HAY fluid from k = 7 to k = 0.5 are presented in Table 
1. Indeed, its values are far from being constant. 



II. CS LAW APPLICATION ANALYSIS 

In order to investigate the applicability of the CS law 
in the case of HAY fluid, we have analyzed the available 
simulation dat a 11 ' 13 ! 14 ] 15 and performed some additional 
Monte Carlo (MC) simulation for relatively long attrac- 
tion tail, k = 1.8 and k = 1.5. Our simulation study has 
been performed in the canonical ensemble as described in 
Refsj n i 13 ' 14 . For the phase diagram definition and sur- 
face tension calculations, MC simulations were performed 
on a parallelepiped cell with dimensions L x — L y = 12, 
and L z > 45. The pressure of the supercritical HAY 



2 



fluid was calculated in a cubic simulation cell, L x = 12. 
In each simulation, the value of cut-off radius r cut of the 
potential (JTJ) was selected to be 5. Such value of r cut as 
well as the system size have been proved to be sufficient 
for obtaining accurate result s 13 ' 14 . Our new results to- 
gether with the previously published results, are given in 
Tables 1 and 2. 

The critical parameters for the HAY fluid were cal- 
culated by using the rectiliniear diameter law and the 
universal value of critical exponent (3 — 0.325. The crit- 
ical pressures were estimated on the base of Clausius- 
Clapeyron equation™-, and together with other critical 
parameters are given in Table 3. 

In fig. 1 we present the dependence of the critical HAY 
fluid density p c and pressure P c , on the inverse critical 
temperature T~ and critical temperature T c , respec- 
tively. Besides our results for the attractive range varying 
from k — 7 to k = 1.5, in the left panel of the figure we 
have used the critical data for k = 0.5, and 1.0 estimated 
from the recent work of Caillol et al^. All the presented 
data indicate a clear linear dependencies described by the 
following equations: 

P c = 0.0228 + 0.0742T C , (3) 
p c = 0.2534 + 0.071-^. (4) 

-L c 

Such liner behavior leads to the constant value of 
the critical compressibility factor, Z c = P c /(p c Tc) = 
0.3 ± 0.01 in the range of k considered (see also Table 
3). This value of Z c coincides surprisingly well with its 
universal value reported in the literature for different real 
substances 1 ^. Actually, in some cases, this value is being 
used to estimate the critical pressure^ 2 -. 

In the left part of the Fig.l we present also the results 
of the so-called self-consistent Ornstein-Zernike approxi- 
mation (SCOZA) reported in Ref A 15 . As seen, there is a 
linear dependence of p c vs 1/T C , too. However, the pre- 
diction of the theory deteriorates at the short attractive 
Yukawa-tails, k > 4. This tendency has been discussed 
in Ref,-, where the authors have addressed some tenta- 
tive reasons of the SCOZA inaccuracy when the range 
of the pair potential narrows. Namely, it was supposed 
that modification of the closure condition equation may 
be needed. 

We have also performed preliminary analysis of the 
available critical data for the fluids with (i) Sutherland- 
like attractive pair interactions, which vary like r~ 3 ~* 
(with t = 3, 1, and 0.1-2S), and (ii) Mie n-6 (7 < n < 32) 
potential 21 . In both cases, there is also a linear depen- 
dence between p c and 1/T C , although it is slightly differ- 
ent from the relation For the case of the Mie fluid, a 
linear function P c = P C (T C ) is presented in the right side 
of Fig.l. It is interesting to note, that there is a region 
where critical parameters of HAY and Mie fluids almost 
coincide. 

On the other hand, the widely used square-well pair 
potential shows, according to SCOZA predictions 2 ^,, com- 



pletely nonlinear behavior of the p c — /9 c (r c _1 ), as can be 
seen in Fig.l. 

If p L and p G denote the liquid and the vapor in mu- 
tual equilibrium at the temperature T, respectively, then 

L G 

according to CS law one should expect p B ] to be univer- 
sal functions of Tr. In Fig. 2 the liquid- vapor coexistence 
curves scaled by the critical values of density and tem- 
perature are presented for different k's. It is seen that 
most of the points lie on or near a single curve. More in- 
terestingly, the rectiliniar diameter line, which is defined 

p L + P G T 

^p^ = l + Al--), (5) 

and formed by the simulation results, is very close 
to the argon diameter line, with A = 3/4, depicted as 
dashed line in Fig. 2. 

Since the surface tension, 7 does not have a reduced 
analogue like Tr or Pr, because jr = 7/7C is not de- 
fined (7 C = 0), the CS law has been successfully ap- 
plied for some substance o 17 ' 20 assuming the universality 
of 7 r = "//(pc^ 3 T c ) with respect to Tr. As we show in 
Fig. 3, all the surface tension simulation data for the HAY 
fluid with k = 7,5, 1.8 and 1.5 map onto the same master 
curve if presented in units of j r . Such CS behavior of the 
reduced surface tension may answer the questions about 
the precision of some recent approximate density func- 
tional theorie s 27 ' 28 ' 29 , which predict unexplained nonli- 
nar trends of j r = 7 r (T^) function. 

The surface tension calculations have been performed 

through the pressure tensor—, therefore, it was natural 

2/3 

to suppose that p r = p/p c (instead of pr) should be 
used in the CS law application for the calculation of the 
reduced pressure, Pr. We have verified such hypothesis 
and the results are presented in Fig. 4 and Table 2. As 
seen, there is a perfect coincidence of Pr for HAY fluid 
with k = 1.8 and 7, at the three reduced temperatures 
considered. 



III. CONCLUSIONS 

This work presents the detailed analysis of the rela- 
tions between the critical parameters of HAY fluid with 
variable range of attraction. Such analysis is based on the 
recently reported accurate simulatio n 11 ' 13 ' 14 ' 15 and inte- 
gral equation data&i 5 -. We show that there is a linear 
relationship between critical pressure and critical tem- 
perature, as well as, critical density and inverse critical 
temperature for HAY fluid for various interaction ranges. 
We have also detected such linearity for the bulk criti- 
cal parameters of n-62i and Sutherland-like 2 ^ potentials 
with varying attractive tail length. We have shown that 
reduced pressure and surface tension of HAY fluid obeys 
the corresponding state law in the range 7 < K < 1.5 
if the re-scaling critical number fluid density p c is re- 
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placed by p c ■ Finally, our findings permit us to con- 
clude that phase diagram, pressure and surface tensions 
of the HAY fluid can be successfully obtained applying 
the properly modified corresponding state theory. Thus, 
we introduce a new criteria for the estimation and im- 
proving of the theoretical approaches^' 7 -' 8 -' 9 - ! 27 ' 28 ' 29 , i.e. 
now the predictive accuracy of different theories can be 
verified by applying of CS law. It is very probable, that 
other fluids interacting via an attractive nonconformal 
potentia l 21 ' 22 ' 25 can be studied with the same CS ap- 



proach. To verify such hypothesis additional accurate 
simulations are desirable. 
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TABLE I: Canonical Monte Carlo data for the phase equilib- 
rium and surface tension of the HAY fluid with n = 1.5 and 
1.8. 



T 


pL 


PG 


7 


K = 1.5 


1.10 


0.730 


0.0207 


0.76020 


1.15 


0.702 


0.0290 


0.633is 


1.20 


0.665 


0.0415 


0.51520 


1.25 


0.633 


0.0525 


0.390i5 


1.30 


0.590 


0.0725 


0.271io 


1.35 


0.545 


0.0930 


0.165i5 


K = 1.8 


0.90 


0.754 


0.0220 


0.612is 


0.95 


0.715 


0.0313 


0.475i5 


1.00 


0.672 


0.0463 


0.360i 


1.05 


0.624 


0.0650 


0.241io 


1.10 


0.565 


0.0880 


0.140io 


1.13 


0.530 


0.1060 


O.O8810 


1.15 


0.500 


0.1432 


0.050 10 



TABLE III: Critical parameters of the HCY fluid with differ- 
ent values of the range parameter k. 



K 


R of 




Pc 




P 




U.u 


15 


D.yoo 


U.ZOU 


-0.4 / 4 






1 n 


15 


a ah 

Z.40U 


U.zoU 


-o.zyu 






1 K 




1.401 


U.oU4 


-O.Z ( ( 


U.lzo 


U.oUU 


1 Q 




l.loU 


U.olo 


£919 
-O.Z1Z 


U.ll^i 


n Qnn 
U.oUU 


l.O 


15 


1 1 GO 


n 017 
U.ol * 


£ 1 99 
-0.1ZZ 






l.O 


11 


1.1/ 1 


n 010 
U.olo 


£ 9/1 9 
-O.Z4Z 






9 m 
z.u 


14 


l.UOU 


u.ozz 


f\ 9^9 

-u.zoz 


U.1UZ 


n Qn^ 

U.oUO 


2.5 


14 


0.840 


0.336 


-6.206 


0.086 


0.303 


3.0 


14 


0.721 


0.356 


-6.028 


0.076 


0.296 


4.0 


13 


0.581 


0.380 


-5.795 


0.066 


0.302 


5.0 


13 


0.500 


0.393 


-5.684 


0.060 


0.305 


6.0 


13 


0.448 


0.412 


-5.590 


0.055 


0.298 


7.0 


13 


0.414 


0.422 


-5.402 


0.053 


0.302 



Figure Captions 



TABLE II: Canonical Monte Carlo results for the supercrit- 
ical pressure, P, of the HAY fluid with k — 7 and 1.8. 





k = 7.0 


K = 1.8 


p 


T = 0.432 


0.50 


0.828 


T = 1.0 


1.425 


2.36 


0.2 


0.052 


0.077 


0.188 


0.130 


0.193 


0.487 


0.4 


0.070 


0.136 


0.478 


0.160 


0.347 


1.280 


0.5 


0.077 






0.210 


0.530 


2.040 


0.6 


0.090 


0.235 


0.900 


0.450 


0.962 


3.420 


0.7 


0.128 


0.351 


1.400 


1.000 


1.855 


5.590 


0.8 


0.240 


0.580 


2.200 


2.400 


3.631 


9.450 


0.9 


0.530 


1.090 


3.500 









Fig.l Critical density p c as a function of inverse critical 
temperature XjT 1 (left panel); critical pressure P c as a 
function of critical temperature T c (right panel). The 
results of the HAY fluids are presented by fill circles 
(MC results), and dotted line (SCOZA^ 15 ). The squares 
depict Molecular Dynamic simulation data of Ref^ for 
n — 6 potential; dashed line represents SCOZA theory 
predictions for the square-well fluid^. 

Fig. 2 Reduced vapor-liquid coexistence curves for 
HAY fluid at different ranges of attractions: n = 7, and 
5 (R e f.i£); k = 1.8 and 1.5 (present work); k = 1 (Ref^). 

Fig. 3 Reduced surface tension 7 r as a function of 
reduced temperature Tr at different ranges of attrac- 
tion. The meaning of the symbols is the same as in Fig. 2. 

Fig. 4 Reduced pressure Pr as a function of reduced 
HAY fluid density p r — p/pV^ at three reduced tem- 
peratures, T R = 1.042, 1.208, and 2.0. Two ranges 
of attraction, k = 7 (squares) and 1.8 (circles) are 
considered. 
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